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AUTHOR'S  SUMMARY 

A  procedure  is  presented  to  calculate  the  compressible  inviscid 

unsteady  transonic  flow  over  an  airfoil,  which  oscillates  sinusoidally  in 

pitch.  In  order  to  treat  precisely  boundary  conditions  on  the  oscillating 

airfoil  surface  and  at  infinity,  the  exterior  of  the  airfoil-shaped  contour 

in  the  physical  plane  is  mapped  onto  a  rectangle  in  a  computational  plane. 

The  two-dimensional  unsteady  Euler  equations  are  solved  there  by  the 

Lax-Wendroff  finite  difference  scheme  with  artificial  viscosity.  Test 

calculations  were  made  for  the  unsteady  flows  over  the  Joukowski  airfoil 

and  the  NACA  0012  airfoil  oscillating  in  pitch,  in  order  to  obtain  several 

individual  flow  patterns.  The  resulting  unsteady  pressure  distributions, 

shock  wave  locations,  etc,  are  presented.  Furthermore,  the  unsteady 

numerical  results  obtained  by  this  procedure  for  the  NLR  7301  airfoil  and 

the  NACA  64A010  airfoil  are  compared  with  the  experimental  ones  by  Tijdeman 

and  Davis,  respectively.  -  - 
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1  INTRODUCTION 

Studies  on  transonic  wings  have  become  active  recently  for  improving  transport 

efficiency,  such  as  in  aircraft  fuel  economy.  This  study  relates  to  unsteady  two- 

dimensional  transonic  airfoils  which  forms  part  of  the  above  studies.  The  transonic  flow 

field  about  an  airfoil  reacts  sensitively  to  a  slight  change  in  boundary  conditions.  For 

instance,  a  change  in  the  angle  of  attack  by  about  1°  results  in  changes  in  shock  wave 

position  and  airfoil  surface  pressure  distribution,  and  great  changes  in  air  force  and 

moment.  It  is  therefore  necessary  to  clarify  such  unsteady  air  force  phenomena  for 

commercialising  a  transonic  airfoil  section.  Examples  of  means  for  this  purpose  are  wind 

tunnel  experiments  and  numerical  experiments.  Examples  of  wind  tunnel  experiments 

currently  carried  out  are  mainly  on  unsteady  boundaries,  such  as  an  airfoil  oscillating 

sinusoidally  in  pitch,  an  airfoil  with  a  flap  oscillating  in  pitch,  and  an  airfoil  with 

time-varying  thickness.  There  are  experimental  measurements  which  have  been  reported  by 
1  a 

Tijdeman  Davis  ,  and  others. 

In  numerical  experiments  too,  a  number  of  methods  of  calculation  for  this  flow  by 
differential  analog  have  been  proposed.  Numerical  experiments  have  been  restricted  to 
inviscid  flow  and  the  transonic  smal 1-disturbance  equation  or  the  full  potential  equation 

or  the  Euler  equations  have  been  taken  as  governing  flow  equations.  Ballhaus  and 

3  ... 

Goorjian  proposed  a  method  of  calculation  using  an  alternating  direction  implicit 

scheme  under  a  low  frequency  approximation,  and  Yu,  Sec-bass  and  Ballhaus'*  proposed  an 
improvement  on  it.  These  methods  can  bo  used  only  where  the  unsteady  disturbance  is  small. 
It  is  impossible  to  obtain  an  accurate  solution  in  the  neighbourhood  of  the  leading  edge 
by  transferring  the  condition  of  tangential  flow  on  the  upper  and  lower  surfaces  to  the 
mean  chord  line  because  of  the  singularity  generated  at  the  nose  of  the  chord  line.  This 
trend  is  particularly  acute  in  the  case  of  a  blunt-nosed  airfoil.  Nevertheless,  because 
of  the  simplicity  in  handling  the  equations  and  boundary  conditions,  these  procedures  are 
inexpensive  and  are  much  quicker  than  other  procedures.  It  is  therefore  possible  to  use 
them  for  calculations  in  a  large  number  of  cases  with  the  use  of  a  mean  angle  of  attack, 
main  stream  Mach  numbers,  and  airfoil  shape,  etc,  as  parameters,  and  to  find  the  trend  of 
variation  in  the  flow  state  corresponding  to  a  variation  in  each  parameter.  In  Japan  too, 
Isogai,  rt  at  developed  a  simple  program  for  a  small-disturbance  potential  equation  with¬ 
out  the  low-frequency  approximation.  With  a  view  to  reducing  inaccuracies  resulting  from 
the  use  of  this  equation,  Isogai^’^  solved  the  full  potential  equation  using  a  semi- 
implicit  time-dependent  difference  method  in  a  cartesian  coordinate  system  and  satisfied 
the  tangential  flow  condition  on  the  fixed  curved  surface  (henceforth  called  the  'mean 
airfoil  surface')  which  coincides  with  the  mean  position  of  the  oscillating  airfoil 

g 

instead  of  the  mean  chord.  Chipman  and  Jameson  carried  out  a  time-dependent  coordinate 
transformation  for  a  sharp-nosed  airfoil,  such  as  a  circular-arc  airfoil,  so  as  to  be 
able  to  give  the  tangential  flow  condition  at  the  true  airfoil  position,  and  solved  a 
full-potential  equation  in  conservation  form.  However,  the  range  of  application  of  these 
potential  equations  is  limited  to  flows  which  contain  only  weak  shock  waves  because  of 
their  isentropic  property.  For  capturing  a  discontinuity  irrespective  of  strength,  it 
is  necessary  to  solve  the  Euler  equations.  The  use  of  these  equations  involves  a  high 
cost  of  computation  and  therefore  not  many  results  have  been  published.  However,  it  is 
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possible  Co  ascertain  the  effectiveness  of  numerical  calculation  by  the  above-referred-to 

two  simpler  equations  by  comparison  of  results  with  those  obtained  from  the  Euler 

9 

equations.  Beam  and  Warming  used  an  implicit  difference  scheme  for  solving  the  Euler 

equations,  but  good  accuracy  could  not  be  expected  for  a  blunt-nosed  airfoil  because  the 

boundary  conditions  were  given  at  a  mean  chord  position.  Magnus  and  Yoshihara^  1-1  gave 

this  condition  at  a  mean  position  on  an  oscillating  airfoil  surface,  and  Lerat  and 

Sides'"’’^  gave  the  same  condition  on  a  true  airfoil  surface  by  a  time-dependent 

coordinate  transformation  and  sought  the  solution  with  the  use  of  an  explicit  difference 

scheme.  For  a  perfect  prediction  of  air  forces,  the  hypothesis  of  inviscidity  should  be 

removed,  in  order  to  enable  the  handling  of  such  viscous  phenomena  as  interaction  between 

1  3 

an  unsteady  shock  wave  and  boundary  layer,  etc.  For  this  purpose,  Magnus  and  Yoshihara 
used  a  method  of  approximating  the  displacement  effect  of  a  boundary  layer,  but  the 
possibility  of  calculations  based  on  the  Navier-Stokes  equations,  which  incorporates  a 
viscous  flow  model,  is  possible.  The  calculation  tried  out  by  Chyu  and  Davis'^  for  by¬ 
passing  the  current  high  cost  of  computation  has  too  coarse  a  grid  spacing.  Satisfactory 
calculations  are  likely  to  be  made  when  the  cost  comes  down. 

This  study  proposes  a  method  of  calculation  for  a  flow  about  an  airfoil  which 

oscillates  sinusoidally  in  pitch,  using  the  unsteady  Euler  equations.  This  method  is  a 

18  19 

variation  of  the  method  of  calculation  for  steady  flow  proposed  earlier  ’  by  the  author, 

by  improving  the  mapping  function  for  grid  construction,  in  order  to  improve  calculation 

efficiency,  and  by  extending  it  to  enable  the  handling  of  unsteady  boundary  conditions. 

It  will  be  assumed  that  the  amplitude  of  pitching  is  small,  and  that  the  boundary  condi- 

1  2 

tions  are  given  at  the  mean  airfoil  surface  position.  Magnus  and  Yoshihara  carried  out 
calculations  in  a  finite  plane  exterior  to  the  mean  airfoil  surface  and  superposed  several 
types  of  Cartesian  grids  varying  in  spacing  on  a  physical  plane  with  a  curvilinear  grid 
near  the  airfoil  surface.  The  difference  scheme  used  was  an  explicit  two-step  finite 

difference  scheme  of  the  Lax-Wendroff  type.  In  this  paper,  the  following  approaches, 

.  .  20 
which  differ  from  the  above,  are  tried  out 

First,  the  infinite  plane  exterior  to  the  mean  airfoil  surface  is  mapped  on  to  a 
rectangle  to  form  a  computational  plane,  in  order  to  facilitate  the  handling  of  airfoil 
surface  boundary  conditions  and  infinity  conditions.  On  the  computational  plane,  only 
one  uniform  grid  is  spanned.  In  the  vicinity  of  a  shock  wave,  a  fine  and  uniform  grid 
is  superposed  when  required.  Numerical  solution  of  the  unsteady  Euler  equations  at  a 
mesh  point  will  be  obtained  from  an  application  of  the  time-split  Lax-Wendroff  scheme, 
with  an  artificial  viscosity  term. 

A  program  will  then  be  made  using  this  method  of  calculation,  and  an  examination 
will  be  made  of  the  possibility  of  obtaining  several  types  of  unsteady  flow  patterns  by 
using  a  Joukowski  airfoil  and  the  NACA  0012  airfoil.  Then,  numerical  solutions  for  the 
NLR  7301  airfoil  and  the  NACA  64A010  airfoil  also  will  be  obtained.  They  will  be 
compared  with  the  results  obtained  by  Tij deman  and  Davis. 
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2  SETTING  I’P  THE  PROBLEM 

The  Euler  equations  of  unsteady  two-dimensional  flow  are  used  as  lundamental 

differential  equations  and  can  be  represented  in  the  following  conservative  form  on  a 

1 8 

physical  plane  x-y.  This  equation  has  previously  been  non-dimensionali sed  in  reference 
by  main  stream  density  and  velocity,  and  chord  length. 
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a  relation  c  =  *'(yp)  /p  holds;  where  t  =  time,  x,  v  rectangular  cartesian  coordinates, 
c  =  density,  u  =  velocity  component  in  x-direction,  v  =  velocity  component  in  y-direction, 
e  =  total  energy  per  unit  volume,  p  =  pressure,  c  =  speed  >f  sound,  y  =  specific  heat 
ratio  of  perfect  gas  (presumed  to  be  1.4). 

The  vector  representation  of  equation  (1)  is 


/  P  \ 


where  W 


W  +  F  (W)  +  G  (W)  =  0 

t  x  y 


The  problem  of  determining  an  inviscid  flow  around  an  unsteady  airfoil  can  be 
replaced  numerically  by  the  problem  of  solving  equation  (2)  subject  to  the  following 
three  boundary  conditions. 


(2) 


(1)  Tangential  flow  condition  -  the  flow  velocity  vector  is  tangent  to  the  unsteady 
airfoil  surface. 


(2)  Condition  of  the  trailing  edge  -  the  Kutta-Joukowski  condition  (the  pressures  on 
the  upper  and  lower  surfaces  at  the  trailing  edge  are  identical)  in  the  absence  of  a 
shock  wave  at  the  trailing  edge  and  the  Rankine-Hugoniot  relation  holds  in  the  presence 
of  a  shock  wave  at  the  trailing  edge. 

(3)  Infinity  conditions  -  at  upstream  infinity,  conditions  are  free  stream  conditions, 

that  is  W  =  Wm  (hereinafter  the  subscript  <*  will  be  used  for  indicating  the  free 

2 

stream),  where  =  1,  *  cos  a^,  v^  =  sin  a^,  p^  =  l/(y(0  where  is  the  mean 

angle  of  attack. 


b 

At  downstream  infinity  and  the  main  tream  subsonic,  shock  wave  can  be  divided  into 
downwash  and  the  rest,  then  p  =  p^,  V  a  ?^,  OU/Sz^)  =  0,  Op/Sz^)  =  0  (z^  *  direction 
of  the  main  stream;  U  =  velocity  spectrum)  hold  for  the  former  and  the  flow  in  the  rest 
is  the  main  stream.  Where  the  main  stream  is  supersonic,  the  downstream  solution  depends 
upon  the  upstream  values,  thus  no  conditions  are  necessary. 

In  this  study,  the  problem  of  unsteadiness  will  be  confined  to  the  flow  about  an 
airfoil  oscillating  sinusoidally  in  pitch  with  small  amplitude,  and  the  conditions  of 
tangential  flow  will  be  applied  at  the  mean  position  of  the  oscillating  airfoil.  The 
method  of  calculation  using  a  computer  for  solving  this  problem  will  be  discussed  in  the 
next  section. 

SOLUTION  OF  FUNDAMENTAL  DIFFERENTIAL  EQUATIONS 
Coordinate  system  and  mesh  for  calculation 

In  order  to  facilitate  the  insertion  of  the  condition  of  tangential  flow  at  the 
airfoil  surface  and  the  condition  at  infinity,  the  mean  airfoil  surface  and  its  exterior 
on  a  physical  plane  (x,y)  will  be  mapped  on  to  a  rectangle  and  its  interior  on  a  computa¬ 
tional  plane  (5,n).  This  mapping  is  carried  out  with  the  following  three  coordinate 
transformations,  which  are  illustrated  in  Fig  I.  The  description  of  mappings  (1)  and  (2) 
is  brief  because  they  have  been  described  in  detail  in  Ref  13  (pp  3,4). 

(1)  A  point  on  an  airfoil  surface  at  the  greatest  distance  from  the  trailing  edge  of  a 
mean  airfoil  surface  is  called  the  leading  edge,  which  is  selected  as  the  origin.  A 
straight  line  connecting  this  with  the  trailing  edge  is  the  x-axis  (a  closed  trailing 
edge  alone  is  considered).  A  mean  airfoil  surface  having  a  trailing  edge  angle  e  and 
its  exterior  on  its  physical  plane  (z  =  x  -  iy)  is  conformally  mapped  on  to  a  unit  circle 
and  its  exterior  on  o  (=  relu) . 

(2)  This  unit  circle  and  its  exterior  are  conformally  mapped  on  to  the  lower  half  of 
w  (=  x'  +  iy')  plane. 

(3)  By  using  the  following  function,  this  lower  half  plane  is  mapped  on  to  a  rectangle 
and  its  interior  on  the  computational  plane  L  =  £  +  in  . 

x’  =  L8(f.2)/(42  -  4)  {U|  £  2) 

y’  =  n/( l  -  n)  (0  s  n  s  I )  . 

In  the  above  equation,  g  is  a  polynomial,  and  its  coefficient  should  be  fixed 
accurately  by  a  method  described  later,  to  obtain  a  highly  efficient  calculation  of  flow 
pattern  which  can  be  visualised  in  each  case  before  starting  the  calculation  of  flow 
patterns . 

In  Fig  I,  the  wavy  line  and  broken  line  show  the  curve  to  which  the  three  straight 
lines  £  «  *b,  n  =  c  correspond  on  each  coordinate  plane. 

With  this  series  of  mappings,  the  mean  airfoil  surface  is  made  to  correspond  to  the 
section  (n  “  0,  -a  f  s  a) ,  upstream  infinity  to  the  upper  side  of  the  rectangle,  and 
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downstream  infinity  to  both  lateral  sides  of  the  rectangle.  Thus,  all  boundary  conditions 
can  be  given  on  the  four  sides  of  the  rectangle.  The  domain  for  computation  is  this 
rectangle  and  its  interior. 

The  Euler  equations  (2)  on  the  physical  plane  are  transformed  into  the  following  on 
the  computational  plane. 


W  +  F.(WK  +  G.(W);  +  F  (W)’  +  G  (W)  •  =  0 

t  X  y  X  r.  y 


13) 


Differential  coefficients  ,  i,  ,  and  ••  ran  be  determined  for  each  mapping 

function,  using  a  chain  rule  (Ref  18,  p.4).  For  solving  equation  (3)  by  a  finite 
difference  method,  first  place  a  uniform  Cartesian  grid  system  having  widths  -V,  and 
on  the  computational  plane  to  make  lattice  lines  coincide  with  the  four  sides  of  the 
rectangle,  and  represent  the  lattice  points  by  symbols  (i,j)  The  grid  image  on  the 
corresponding  physical  plane  is  also  orthogonal.  As  will  be  shown  later  in  examples,  the 
j  group  of  lattice  lines  are  higher  in  density  as  they  approach  the  airfoil  surface, 
whereas  the  density  distribution  among  the  i  group  of  lattice  lines  is  designed  to  be 
adjustable  by  means  of  the  function  g  in  the  mapping  in  equation  (3).  Taking  into 
cons ide rat  ion  the  pattern  of  a  predicted  flow,  g  is  fixed  to  make  the  grid  image  density 
on  the  physical  plane  higher  in  an  area  where  flow  variation  is  great  (in  the  vicinity  of 
shock  wave,  leading  edge  and  trailing  edge).  However,  arbitrary  points  j  and  j  on 
the  :',-axis  downstream  of  the  trailing  edge  must  correspond  to  the  same  point.  The 
function  g  is  made  a  function  of  ,  Because  of  this,  the  lattice  distributions  on 
the  upper  and  lower  surfaces  of  an  airfoil  are  similar  in  shape.  However,  where  the 
flow  patterns  on  the  upper  and  lower  surfaces  differ  considerably,  for  instance,  in  the 
presence  of  a  shock  wave  only  on  the  upper  surface,  the  grid  in  its  vicinity  must  be  made 
fine,  but  there  is  no  need  to  do  this  on  the  lower  surface.  For  maintaining  computational 
efficiency  in  such  a  case,  a  grid  (width  was  formed  on  first  fixing  g  to  match  the 
flow  pattern  on  the  lower  surface,  and  a  finer  grid  (width  .':./4)  was  superposed  on  that 
part  of  the  grid  system  corresponding  to  the  vicinity  of  the  shock  wave  on  the  upper 
surface.  An  automatic  superposition  of  that  section  with  the  l/4-width  grid  was  possible 
by  the  inclusion  of  a  routine  lor  designating  (i  ,i0).  With  the  selection  of  a  suitable 
g  ,  and  the  superposition  of  the  I /4-width  grid,  an  improvement  in  computational 
efficiency  and  an  accurate  numerical  grasping  of  flow  pattern  variation  were  made 
possible . 

Because  the  boundary  conditions  of  an  oscillating  airfoil  are  applied  at  the  mean 
position  of  the  airfoil  surface  a  time-dependent  grid  system  is  not  required. 


The  definition  of  g:  The  correspondence  between  z  =  x  -  iy  and  w  =  x'  +  iy *  is 
obtainable  from  mappings  (1)  and  (2),  therefore  we  will  try  to  obtain  a  curve  (inclusive 
of  a  mean  airfoil  surface)  corresponding  to  the  x'-axis.  Choose  a  suitable  number  of 


points  (x  .y^),  where  t  =  I,  2,  3 .  L  ,  on  that  curve  and  arrange  the  desired 

lattice-point  subscript  i  for  each  of  these,  while  taking  into  consideration  the 


predicted  flow  pattern.  Then,  find 
of  data  [q,  ,  which 


x'  corresponding  to  (x^.y^)  and 

can  be  calculated  from  x'  and 

s 
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with  a  number  of  least  squares  polynomials  differing  in  degree.  Draw  a  grid  image  lor 
each  degree  on  the  physical  plane,  select  one  which  is  likely  to  be  suitable  lor  flow 
calculation,  and  define  g  as  follows. 


M 


m=0 


.  18 

3 . J  Finite  difference  method  at  an  interior  point 

The  solution  by  finite  differences  for  W  at  the  point  (i,i),  alter  >.'  time  steps 
is  represented  by  .  .  For  mesh  points  in  the  rectangle,  the  following  tiae-spli: 

difference  approximations  are  applied  to  equation  1 3)  <. F i ^  IT). 


WN+!  =  L  L.WN  .  ,  WN  +  ‘  =  V  V.L.l  WN+1 

r,J  -  A  i.j  i.J  '  •  •  r  i.J 


The  differential  operator  L.  in  the  above  equations  is  a  two-step  Lax-Wendrof f  operator 
of  second  order  accuracy. 


L.W.  .  =  W.  .  -  ~  | ( F .  +  (C.  ,  .  -  G.  ,  .):  | 

'  t,J  i.J  i  +  :»J  f-i.J  :  1+ :  >  J  y  [  : 


W.  ,  .  =  *(W.  .  .  +  W.  .)  -  I  (F ,  .  -  F .  .)  i  +  (G.  .  •  “  G.  .  )  < 

i+i.J  ‘  i  +  l,J  i.J  2 At,  j  l+l, j  i.J  xi+ij  1  +  1>J  i.J  >';  +  }, j! 


can  be  defined  similarly.  The  von  Neumann  condition  for  L.L  L  L.  can  be  given 
by  the  following  equation 


1  8 


At  <  min 


(  &  ■'  - _ ) 

V  I  u5  +  Vt,  :  +  c  ’  iu-1  +  V”  +  cl 
\  x  y  -  x  v 1  ' 


=  AT  . 


For  obtaining  a  value  for  the  magnitude  of  time-step  for  use  in  actual  calculations , 
multiply  AT  by  a  suitable  scalar  (less  than  1),  so  that  the  stability  criterion 

is  satisfied.  The  artificial  viscosity  operator  V.  can  be  defined  as  follows: 


V  W.  . 
'■  i.J 


w.  .  +  \r  •Jlu1:  ,  .  -  u'r  .  1  <w .  .  .  -  w.  .)  -  !ut:  .  -  if:  .  .  •  (w.  .  -  w.  .  .)}  , 

i,j  i)  1+1, j  i.j  i+i,j  i.j  i.j  i-i, j  i,j  l-i, j  | 


+  vC 
x  y 


is  defined  similarly.  The  quantities  A  and  \  are  suitably  chosen  to  ensure 
the  numerical  stability  of  the  solution  with  artificial  viscosity  terms. 


L 
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3.3  Method  of  handling  boundary  conditions 

As  was  said  in  section  3.1,  the  whole  of  the  airfoil  surface,  trailing  edge,  and 
infinities  correspond  to  the  four  sides  of  the  rectangle  on  the  computational  plane. 

Thus  a  method  for  obtaining  numerical  solutions  which  satisfy  the  boundary  conditions 
discussed  in  section  2  will  be  described. 

(1)  Two  sections  (a  <  2  and  -a  -  £,  '■  -2)  on  the  base  (n  =  0)  :  They  correspond 

Co  the  same  curve,  which  is  not  a  boundary,  on  the  physical  plane  in  Fig  I  (the  wavy- line 
sections  in  Fig  II).  Thus  the  method  relating  to  the  internal  point  discussed  in 
section  3.2  will  be  applied  with  the  relation  W((.)  =  W(-f.)  assumed  to  hold. 

(2)  Section  i;  '  a  on  the  base  (n  =0):  This  corresponds  to  mean  airfoil  surface. 
The  condition  of  sinusoidal  oscillation  of  the  solid  airfoil  surface  must  be  satisfied. 

In  other  words,  the  velocity  component  of  a  flow  normal  to  the  airfoil  surface  must  equal 
the  component  of  velocity  of  the  surface  in  the  same  direction.  V.’  on  the. airfoil 
surface  in  this  calculation  is  obtainable  from  the  following  representation  (j  =  1 
indicates  the  base) . 


0L  L.W1?" 
h  9  1,1 


W* 


+  2  =  OV.L.iu1?*! 

c  n  1,1 


*>  *  A  A 

In  the  above  equation,  l  can  be  defined  as  l  U.  =  h(L  W.  .,  L  W.  .)  where  h  is 

^  1)1  n  i  ^  i  j  3 

an  extrapolation  function  relating  to  the  distance  between  corresponding  points  on  the 
physical  plane  of  grid  points  j  =  (,  2,  3.  First,  the  intermediate  quantity  W  will 
be  obtained  from  the  following  equation: 


W?+?  =  VfLi^! 

1,1  U  u,l 


The  velocity  components  in  W  are  replaced  by  the  tangential  component,  v  and 
the  normal  component  vnor  at  that  time.  This  quantity  W  can  be  updated1"*  to  W  by 
an  oscillation  operator  0  which  is  used  in  many  operations  to  be  discussed  later. 

We  will  now  consider  (Fig  III)  a  coordinate  system  X  -  Y  (X  =  direction  of  the 

airfoil)  which  is  fixed  to  the  moving  airfoil,  in  addition  to  coordinate  system  x  -  y  , 

on  the  physical  plane.  When  the  airfoil  is  in  its  mean  attitude,  they  coincide  with  each 

other.  Let  us  assume  that  an  airfoil  is  oscillating  sinusoidally  in  pitch  about  a 

central  axis  of  coordinates  (x  ,y  )  =  (X  ,Y  ),  with 

c  c  c  c 

a(t)  =  “q  +  8(t)  ,  8(t)  =  Aa  sin  2kt 

where  a(t)  =  angle  of  attack,  =  mean  angle  of  attack,  6(t)  =  angle  of  pitching, 

Aa  «  amplitude  of  pitching,  k  =  frequency  parameter. 

The  frequency  parameter  k  is  defined  by 


Ill 


k  = 

where  i  =  airfoil  chord  length,  -  '  =  angular  frequency,  and  1  =  specif  ol  main  stream. 
The  po.ltion  of  point  m  (coordinates  corresponding  to  a  gnu  tixed  to 


ving  airfoil  in  the  : 

c-y  system  at 

the  time  t  is 

ob  Lai  nab  1 e 

from 

t  hi*  uPjUd  t  i  or.  s 

x  U) 

=  X  *  (X 

-  X  )  COS  :■  1  t  ) 
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u  > 

'  m  erne  no 

The  velocity  component  in  the  outward  normal  direction  at  pe  nt  m  on  the  a i r ; c i  1 

is  obtainable  from  the  velocitv  components  lx  ,v  )  in  the  Cartesian  system  bv  means  of 

m  "  m 

the  equation 

v  =  -  x  sinfc  -  :)  +  v  eos(’  -  .) 

nor  m  n  ■  m  m 

where  v  is  the  angle  between  the  tangent  to  the  airfoil  surface  and  the  airfoil  chord 

measured  positive  in  the  clockwise  direction.  The  velocitv  component  v  in  a  direc- 

nor 

tion  normal  to  the  airfoil  surface  at  the  time  t  ,  which  has  been  obtained  earlier,  is 

adjusted  to  agree  with  the  velocitv  component  v  of  a  plane  isentropic  wave  which 

nor 

propagates  in  the  normal  direction  on  the  outside  of  the  airfoil.  This  adjustment  causes 

variation  in  density  and  pressure,  but  it  does  not  change  the  tangential  component  of 

velocity  v  .  The  adjusted  density,  pressure  and  flow  velocitv  (v  ,v  1  are  used 
tan  •  *  -tan  nor 

for  updating  W  . 

'31  Two  points  (;,  =  "a,  r  =  0)  on  the  base:  These  points  correspond  to  the  same  point, 
the  trailing  edge,  on  the  physical  plane.  Although  the  conditions  differ  according  to 
the  presence  or  absence  of  a  shock  wave  at  the  trailing  edge,  as  was  said  earli.r,  a 
J i scon t i nuous  surface  is  handled  as  a  continuous  surface  with  a  steep  slope  in  the  finite 
difference  method.  Thus,  for  convenience,  in  the  presence  of  a  mesh  poin1-  at  the  trail¬ 
ing  edge,  numerical  solutions  for  the  upper  and  lower  surfaces  are  obtained  and  these  are 
then  averaged.  In  the  absence  of  such  a  mesh  point,  use  the  method  relating  to  the 

internal  point  discussed  in  the  previous  section  at  two  mesh  points  (L  =  , ,  -  =  0) 

nearest  to  the  base  of  the  rectangle  and  average  their  numerical  solutions. 

(.4)  <2,  n  =  lj  on  the  top  side:  This  corresponds  to  upstream  at  infinity, 

thus  it  always  gives  main  stream  . 

(5)  ^ i  £.  r2,  0  ••  n  <  lj  on  two  lateral  sides:  This  corresponds  to  uownstream  at 

infinity.  Where  the  main  stream  is  supersonic,  the  numerical  solution  at  a  neighbouring 
grid  point  within  the  rectangle  is  used  at  the  boundary  point.  In  other  words,  extrapola¬ 
tion  with  zero-order  precision  is  carried  out.  On  the  other  hand,  if  the  main  stream  is 
subsonic,  the  condition  at  infinity  discussed  in  section  (3)  should  be  imposed.  However 
it  is  difficult  to  use  that  condition  in  the  present  calculation,  so  the  following 
numerical  approximation  is  made.  Whether  or  not  a  grid  point  (IE,j)  at  C  =  :2  is 
downstream  of  a  shock  wave  is  determined  (Fig  IV)  bv  whether  the  value  of  entropv  S 
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ROUT 1  Nr.  OF  NUMERICAL  vALCULAl ION 

Program  'MESH'  for  terming  a  mesh  with  trie  selulioi.s  cisi  ussec  it.  ::>e  pr,  v: . 
chapter,  program  'PLOW'  for  finding  a  flow  about  an  airfoil  by  the  tir.ite  di:  ten > 
method,  and  program  'RESULT'  processing  results  based  or.  V  ,  with  graphs,  etc,  were 
prepared.  These  programs  were  written  in  the  code  for  a  FACOM  2 '0-75  computer  with  art 
array  processor  unit. 

The  calculation  of  the  case  of  an  airfoil  in  :  main  stream  ot  given  Mach  number  M 

oscillating  in  pitch  with  angle  of  atlack  ( t  '  =  i  +  i  sin  2k  t  about  the  central  axis 

IX  ,Y  1  was  carried  out  in  the  following  sequence  of  steps, 
c  0 

(. I )  Preparation  of  the  mesh:  The  function  g  was  determined  with  program.  Mid SH  by  : 

procedure  discussed  in  section  3.  1  .  It  was  decided  whether  to  superpose  the 

grid  and  if  so  the  position  at  which  the  superposit ion  was  to  be  made  was  determined. 

(2)  Calculation  of  a  steady  flow  which  becomes  a  starting  solution  in  the  calculation 
of  unsteady  flow:  An  unsteady  flow  about  an  airfoil  placed  at  an  angle  of  attack  i„  it. 
a  main  stream  of  Mach  number  M  was  found  by  the  following  procedure.  Seek  asymptoti¬ 
cally  steady  solution  by  applying  the  program  FLOW  on  a  suitable  starting  value  with 
amplitude  d.i  =  0  .  (.Described  in  detail  in  Ref  18.) 

(J)  Calculation  of  a  quasi-steadv  flow:  Although  net  always  necessary,  it  is  advisable 
to  have,  before  starting  the  calculation  for  an  unsteady  flow  (k  *  01 ,  a  calculation  f.r 
the  corresponding  quasi-steadv  flow  (k  =  0,  oscillation  is  infinitely  slow)  tor  refertti.e. 
That  is,  asymptotical lv  steady  solutions  W  are  calculated  for  both  maximum  angle  of 
attack  (Uq  +  /.a)  and  minimum  angle  of  attack  (  < )  of  the  oscillatioti,  as  in  (21 

above.  Convergence  is  obtainable  more  rapidly  from  using  a  starting  solution  obtained  by 
the  method  discussed  in  (2).  It  is  necessary  at  this  time  point  to  examine  whether  a 
mesh  prepared  in  (1)  is  suitable  for  calculations  in  12)  and  (3).  If  unsuitable,  repeat 
the  preparation  in  accordance  with  (1).  If  the  mesh  has  been  found  suitable  for  calcula¬ 
tions  in  (2)  and  (3),  the  calculations  of  the  unsteady  flow  having  an  angle  of  atlack 
a(t)  continually  varying  between  the  above  two  limits  can  probably  be  made 
satisfactorily. 

(4)  Processing  the  results  of  the  quasi-steady  flow  calculation:  Derive  the  following 
from  the  results  of  calculations  in  (2)  and  (3),  using  RESULT. 
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(a)  Prepare  isobaric  charts  with  sonic  lines  tor  three  steady  tlows  '  angle: 


attack 


0 


ja ,  a 


O’ 


-  A  a )  and  diagrams  |C  '  2-.p-  1/(.M“)  |  of  variation  in  airfoil 


surface  pressure  coefficient  Cp  with  chord  length.  Integrate  Cp  to  obtain  the 
coefficient  of  drag  ,  the  coefficient  ot  lift  C^  ,  and  the  moment  coefficient  C^ 

V  diagram  of 
with  chord  length,  where 


(b)  A  diagram  of  variation  in  quasi-steady  pressure  coefficient  AC  (.or 


S  ' 


AC 


E  (cp(a0  -  Aa)  -  LpU0  -  ->]/'-■■• 


ll 


or,  in  the  case  of  a  symmetrical  airfoil,  where  a^  =  0  : 


(5)  Calculations  for  unsteady  flow:  With  the  flow  tor  the  mean  angle  of  attack  obtained 
in  (2)  above,  taken  as  starting  solution  at  t  =  0  ,  an  unsteady  flow  about  an  airfoil 
oscillating  sinusoidally  in  pitch  will  be  obtained  with  the  use  of  FLOW.  Continue  calcu¬ 
lations  until  periodic  solutions  are  obtained  (about  four  periods  are  necessary  for 
coinciding  isobaric  charts  for  all  periods).  During  that  time,  solutions  W  obtained 
for  equidistant  phase  angles  (every  10°  here)  were  stored  in  a  disc  pack  to  be  used  later 
for  processing  the  following  results. 

(6)  Processing  of  results  pertaining  to  unsteady  flow:  RESULT  was  used  for  deriving 
the  following  from  the  results  W  stored  in  (5)  above. 

(a)  Isobaric  charts  with  sonic  lines  at  a  series  of  specified  phases,  and  diagrams 
of  variation  in  airfoil  surface  pressure  coefficient  with  chord  length. 

(b)  Unsteady  pressure  distribution  diagram  (depiction  of  first  harmonic  variation 
of  C  ).  That  is,  to  obtain  CD-  at  a  lattice  point  on  an  airfoil  surface  in  everv 

p  r  i 

r  o 

equidistant  phase  angle  (10  )  from  the  solution  W  for  the  final  period  stored  in  (5), 
and  then  approximate  the  data  (37  in  number)  for  the  period  under  consideration  by  a 
spline  interpolation  and  a  fast  Fourier  transform  to  get  the  equation 

C  (t)  "  Ba  +  AafRefc  ^  sin  2kt  +  Im/c  \  cos  2kt; 

Pi  °i  I  V  Pi>  '  PiJ  \ 


=  B„  +  Aa  I  AC  |  sin(2kt  +  0.) 
°i  '  Pi'  1 


In  the  above  equation,  Re(Cp)  is  the  real  part  of  unsteady  pressure  distribution 

which  is  of  the  same  phase  as  the  movement  of  airfoil,  and  Im(C  )  is  the  imaginary  part 

P 

of  the  unsteady  pressure  distribution  which  is  in  a  quadrature  with  the  motion.  Further, 
| AC  |  is  the  magnitude  of  the  first  harmonic  and  b  is  phase  relative  to  the  airfoil 
motion.  Variation  of  these  four  quantities  with  chord  length  is  shown  in  a  figure. 


(c)  Air  force  coefficient  and  Fourier  analysis  of  the  position  of  the  shock  wave. 

The  position  of  the  shock  wave  (X  of  minimum  C  )  was  calculated  from  C  obtained  during 

P  P 

the  calculation  in  (b)  above.  Tangential  force  coefficient  (C  ) ,  normal  force 

A 


coefficient  (C^,)  ,  pitching  moment  coefficient  ICX^.V^)  centre  and  positive  nose-op  i 

were  calculated  for  each  phase  angle  (10  )  by  integration.  Then,  each  of  these  was 

represented  by  f  and  the  data  (37  in  number)  thus  obtained  for  one  period  were 

approximated  by  the  following  equation,  obtained  by  spline  interpolation  and  fast 

transform.  Numerical  values  of  E„,  A  and  ;  were  obtained 

O  n  n 

N' 

f(t)  =  Bq  +  'a 

n=  1 


D 


A  sin(2nkt  +  ;  )?  . 
n  n  ) 


If  the  magnitude  of  the  second  and  subsequent  harmonics,  A^  ,  where  n  =  2,  i,  ..., 
is  smaller  than  that  of  the  first  harmonic  A  ,  f  changes  sinusoidally  with  respect 
time.  In  other  words,  it  can  be  said  that  a  linear  relation  exists  between  the  displace¬ 
ment  of  an  airfoil  and  its  f  .  Where  this  is  not  so,  these  coefficients  serve  as  an 
indicator  of  the  degree  of  deviation  from  the  sinusoidal  form. 

5  EXAMPLES  OF  NUMERICAL  CALCULATION 

Let  us  illustrate  the  numerical  capturing  of  various  types  of  unsteady  shock  wave- 
phenomena  at  various  main  stream  Mach  numbers,  by  the  calculation  sequence  given  in  steps 
1  to  5.  Experimental  conditions  differ  slightly  fr  m  the  conditions  of  calculation  usee 
here  because  of  the  effects  of  interference  by  wind  tunnel  walls  and  viscosity,  but  cases 

6  and  7  illustrate  the  reproduction  by  calculation  of  an  analogous  flow  pattern  which 
enables  a  comparison  with  the  experiments  of  Tijdeman  to  be  made.  In  case  8,  the 
present  method  of  calculation  are  compared  with  the  results  of  Davis’  experiments,  anc 
of  other  methods  of  calculation. 

The  types  of  airfoil,  the  values  of  parameters  (main  stream  Mach  number,  mean  angle 
of  attack,  amplitude  of  pitching,  reduced  frequency,  position  of  central  axis  of  oscilla¬ 
tion,  number  of  mesh  points,  coefficient  of  time-step  and  coefficient  of  artificial 
viscosity),  the  presence  or  absence  of  superposition  with  AC/4-wide  mesh,  and  numbers 
of  mesh  images  used  in  the  calculation  in  each  are  shown  collectively  in  Table  1 .  The 
numbers  of  the  figures  showing  the  results  of  the  principal  calculations  in  all  cases  ari 
given  in  Table  2.  The  symbols  (  )  indicate  the  numbers  of  the  figures  showing  comparable 
results  obtained  by  experiments  or  other  methods  of  calculation.  All  figures  pertain  to 
a  period  a  number  of  cycles  2kt/2ir  from  t  =  0  have  occurred.  The  characteristics  and 
results  of  the  calculations  will  now  be  examined  case  by  case. 

5 . I  Tests  on  capturing  numerical  values  pertaining  to  various  types  of  unsteady- 

flow  fields 

(Case  1)  Numerical  values  pertaining  to  unsteady  phenomena  involving  the  generation, 
growth  and  annihilation  of  a  weak  shock  wave 

At  the  free  stream  Mach  number  of  0.73,  a  steady  flow  produced  by  the  symmetrical 
airfoil  NACA  0012  does  not  contain  a  shock  wave  when  the  angle  of  attack  is  0°,  as  is 
shown  in  Fig  l-2f,  but  it  does  when  the  angle  of  attack  is  2°.  Under  the  conditions  of 
a  mean  angle  of  attack  0°  and  amplitude  of  pitching  2°,  calculations  at  k  =  0.1  were 
made.  The  development  of  the  flow  was  followed  numerically,  through  the  generation  of  a 


Values  of  parameters  and  of  grids  used  in  calculations 
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comparing  the  present  method  with  experiment  or  other  ea 1 eu 1  at  ion  method. 
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shock  wave,  its  growth  and  decay  on  the  upper  surlace,  followed  by  the  same  events  on  the 
lower  surface  in  every  period.  Periodic  solutions  were  obtained,  as  is  seen  in  Figs  1-1! 
and  1-1  but  they  do  have  a  phase  lag.  It  should  therefore  be  possible  to  superpose  a 
diagram  of  a  flow  field  with  a  random  cycle  diagram  ot  another  flow  field  dillering  by 
0.5  cycle  from  the  first.  This  has  been  proved  (Fig  1-2).  As  will  be  seen  from  a  diagram 
(Fig  1--,)  of  quasi-steady  pressure  distribution  and  diagrams  (Figs  1-5  and  1-61  of 
unsteady  pressure  distribution,  the  absolute  value  of  unsteady  pressure  distribution 
coefficient  is  greater  on  the  forward  part  of  the  airfoil.  This  is  the  part  on  which 
generation,  growth  and  decay  of  shock  waves  take  place. 

lease  2)  Numerical  capturing  of  unsteady  phenomena  in  displacement  of  a  strong 

shock  wave 

Under  the  condition  of  a  free  stream  Mach  number  of  0.72,  a  strong  and  persistent 
shock  wave  is  produced  on  the  upper  surface  of  a  Joukowski  airfoil  at  an  angle  of  attack 
of  0°,  2°  or  A0.  Unsteady  calculations  with  k  =  0. 1  and  0.25  were  made  under  conditions 
of  a  mean  angle  of  attack  2°  and  an  amplitude  of  oscillation  of  2°.  With  either  value  of 
k  ,  the  maximum  amplitude  of  displacement  (Fig  2-5b&c)  of  in  a  period  was  smaller 

than  an  absolute  difference  (Fig  2-5a)  in  between  two  steady  flows  having  angles  of 

attack  equal  to  0°  minimum  and  4°  maximum.  It  will  be  seen  from  Figs  2-3  and  2-4  that  the 
supersonic  region  continues  to  enlarge  even  after  passing  the  time  (planes  C  in  these 
figures)  of  a  maximum  angle  4°  of  attack,  and  continues  to  shrink  after  passing  the  time 
(planes  G  in  the  figures)  of  minimum  0°  angle  of  attack,  indicating  a  phase  lag.  In 
the  unsteady  pressure  distribution,  Fig  2-7  shows  a  pressure  peak  originating  from  the 
displacing  shock  wave.  With  an  increase  in  frequency  parameter  k  ,  the  pressure  peak 
deviates  gradually  from  the  real  part  to  the  imaginary  part.  In  other  words,  with  an 
increase  in  frequency,  periodic  shock  wave  motion  increases  the  time  lag.  It  will  be  seen 
from  Fig  2-8,  which  illustrates  unsteady  pressure  distribution  with  absolute  values  and 
phase  angle,  the  width  and  height  of  a  pressure  peak,  which  is  coupled  to  shock  wave, 
decreases  with  an  increase  in  frequency.  This  was  caused  by  a  reduction  in  amplitude  of 
displacement  with  an  increase  in  frequency. 

Coefficients  up  to  the  third  harmonic  of  tangential  force  (C^),  normal  force  (C^,), 

pitching  moment  (C  )  (about  centre  of  half  chord  and  positive  nose-up),  and  the  position 
M 

of  shock  wave  (S^)  (position  of  lowest  pressure  on  the  upper  surface),  which  were  given 
by  Fourier  analysis  in  section  (6c)  of  the  previous  section,  are  shown  in  Table  3.  At 
either  value  of  k  ,  the  ratio  of  the  second  harmonic  to  the  first  harmonic  in  the 
coefficient  of  normal  force  or  coefficient  of  tangential  force  is  smaller  than  0.06.  It 
is  therefore  indicated  that  the  loci  of  these  coefficients  are  close  to  sine  curves 
incorporating  different  types  of  time  lag,  despite  the  presence  of  an  oscillating  shock 
wave.  However  for  the  coefficient  of  tangential  force  and  the  position  of  shock  wave, 
the  ratios  are  greater  than  0.1.  In  other  words,  it  is  clear  that  these  two  coefficients 
deviate  from  sinusoidal  loci  and  behave  nonlinearly.  With  an  increase  in  the  value  of  k 
from  0.1  to  0.25,  considerable  reduction  in  absolute  values  of  the  first  harmonics  and 
increases  in  their  phase  angles,  except  that  of  ,  were  observed. 


Table  J 
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Fourier  coefficients  for  coefficient  of  chord  force  (C  )  ,  coefficient  of  normal  foi  -e  (C^,)  , 

~  coefficient  of  pitching  moment  (C  )  (about  centre  of  half-chord,  positive  nose-up) 
_  N _ _ _ 

and  the  position  of  shock  wave  (S^) 


(a)  k  =  0. 1 


Mean 

Amplitudes 

Phase  angles 

f 

Bo 

A1 

A2 

A  3 

'  1 

'  3 

cx 

0.0047 

0.0064 

0.0034 

0.0004 

-32.97 

48.77 

116.79 

CY 

1 .0685 

0.2176 

0.0052 

0.0012 

-  1  8 . 34 

17.68 

-156.84 

CM 

0.0813 

0.0292 

0.0016 

0.0004 

-22.02 

0.44 

40.7  1 

sx 

0.6058 

0.0340 

0.0040 

0.0002 

-36.31 

21.94 

17  5.13 

(b)  k 

=  0.25 

Mean 

Amplitudes 

Phase  angles 

f 

Bo 

A. 

A2 

A3 

f  1 

**  o 

;  3 

cx 

0.0024 

0.0040 

0.0026 

0.0004 

-66.79 

29.33 

94.10 

CY 

1.0505 

0.1586 

0.0014 

0.0006 

-15.96 

-34.68 

-145.20 

S, 

0.0858 

0.0270 

0.0002 

0.0000 

-27.61 

33.16 

59.90 

sv 

0.5968 

0.0240 

0.0032 

0.0006 

-70.79 

-60.27 

-149.50 

X 

- 

— 

(Case  3)  Capturing  of  an  unsteady  flow  containing  strong  diagonal  shock  wave  always 
present  at  the  trailing  edge  on  the  upper  surface  (downstream  is  subsonic), 
and  shock  wave  periodically  displacing  downstream  of  the  trailing  edge 
on  the  lower  surface 

At  a  main  stream  Mach  number  of  0.93,  a  steady  shock  wave  is  present  at  the  trailing 
edge  on  the  upper  surface  of  a  Joukowski  airfoil  where  the  angle  of  attack  is  -3°,  0°  or 
3°,  as  shown  in  Figs  3-2  and  3-Aa.  On  the  lower  surface,  the  position  of  another  shock 
wave  moves  a  large  distance  upstream  from  the  trailing  edge  with  an  increase  in  angle  of 
attack.  Where  the  mean  angle  of  attack  is  0°,  the  amplitude  of  pitching  is  3°  and  the 
frequency  is  0.1,  the  unsteady  shock  wave  (Fig  3-3)  on  the  low. r  surface  moves  periodically 
within  a  range  surrounded  by  two  steady  shock  waves  (Fig  3-2)  produced  at  an  angle  of 
attack  i3°.  It  will  be  seen  from  Fig  3-4  that  the  variation  of  C^  within  the  period 
relevant  to  k  =  0.1  is  smaller  than  that  of  of  a  quasi-steady  flow.  As  will  be 

seen  from  Figs  3-5,  3-6  and  3-7,  the  pressure  coefficient  for  unsteady  flow  is  a  flat 
curve  close  to  0  on  the  upper  surface  because  of  an  almost  fixed  shock  wave,  but  the 
absolute  values  vary  widely  on  the  rear  half-chord  because  of  the  large  movement  of  the 
shock  wave. 


18 


(.Ease  ->)  Capturing  numerical  values  of  'fish  tail'  shock  wave  system  and  the  process 
of  its  degradation 

The  appearance  of  the  so-called  'fish  tail'  shock  wave  system  is  well  known  when  the 
free  stream  Mach  number  is  close  to  1.  That  is,  weak  oblique  shock  waves  (flow  changing 
: too  supersonic  to  supersonic)  attach  to  the  trailing  edge  on  the  upper  and  lower  surfaces 
and  merge  with  a  normal  shock  wave  downstream  01  the  trailing  edge.  In  this  triangular 
legion  formed  by  these  three  shock  waves,  there  is  no  great  variation  in  Mach  number.  At 
a  l  ree  stream  Mach  number  of  0.93,  a  shock  wave  system  having  a  similar  structure  was 
obtained  with  a  Joukowski  airfoil  when  the  angle  of  attack  was  -3°.  However,  when  the 
angle  of  attack  was  0°,  a  fairly  small  subsonic  region  appeared  behind  a  steady  shock  wave 
;  •;  the  vicinity  of  the  trailing  edge  on  the  lower  surface.  When  the  angle  of  attack  was 
>  ,  this  subsonic  region  was  much  bigger.  (Fig  4-1,  a  slip  flow  was  present  at  a  position 
where  a  sonic  line  of  a  normal  shock  wave  changed  rapidly  downstream  of  the  trailing  edge) 
An  unsteady  solution  was  then  sought  under  conditions  of  mean  angle  of  attack  0°,  amplitud 
of  pitching  3°  and  frequency  parameter  0,1.  During  the  process  of  obtaining  the  periodic 
solution  shown  in  Fig  4-3,  a  'fish  tail'  shock  wave  system  was  maintained  at  time  points 
G,  H  and  A  (angles  of  attack,  -3°  -*■  -1.5°  -»  0°>,  a  supersonic  region  spread  out  in  the 
vicinity  of  the  trailing  edge  on  the  lower  surface  with  an  advance  in  time  points  to  B, 
i  and  1)  (.1.5°  -  3°  ->•  1.5°),  it  contracted  with  an  advance  to  D,  E  and  F  (1.5°  0°  v 

-1.3°),  and  disappeared  at  time  point  G  (-3°).  As  compared  with  the  previous  case 
(Fig  3-3),  the  shock  wave  on  the  upper  surface  (Fig  4-2)  leaned  considerably  towards  the 
direction  of  Che  chord  line.  The  diagrams  (Figs  4-3,  4-4,  4-5  and  4-6)  of  steady  and 
utv  teadv  pressure  distribution  on  the  airfoil  surface  differed  considerably  from  the 
previous  case  in  the  patterns  of  flow  in  the  vicinity  of  and  downstream  of  the  trailing 
edge,  but  were  the  same  in  trend.  However,  the  peak  in  an  absolute  value  of  unsteady 
pressure  coefficient  caused  by  the  shock  wave  on  the  lower  surface  was  smaller. 

(Ease  5)  Capturing  numerical  values  of  an  unsteady  flow  about  an  airfoil  where  the 
main  stream  Mach  number  is  greater  than  I 

When  the  free  stream  Mach  number  exceeds  1  an  isolated  shock  wave  forms  upstream 
of  a  circular  leading  edge,  which  is  almost  normal  to  the  air  flow  in  the  vicinity  of  the 
chord  line,  thus  it  becomes  subsonic  downstream  of  it.  However,  it  is  well  known  that 
this  is  accelerated  again  to  supersonic  speed  on  the  airfoil  surface,  and  a  weak  oblique 
shock  wave  (its  flow  changing  from  supersonic  to  supersonic)  forms  at  the  trailing  edge. 
Here  again,  in  respect  of  a  Joukowski  airfoil,  similar  steady  flow  patterns  were  obtained 
under  conditions  of  a  free  stream  Mach  number  of  1.4  and  angles  of  attack  of  -5°,  0°  and 
5°  (Figs  5-1  and  5-3a) .  The  supersonic  region  in  the  vicinity  of  the  leading  edge  was 
the  greatest  where  the  angle  of  attack  was  -5°.  In  other  words,  an  isolated  shock  wave 
was  the  smallest.  According  to  the  results  (Figs  5-2  and  5-3b)  of  calculating  an  unsteady 
flow  under  conditions  of  a  mean  angle  of  attack  0°,  amplitude  of  pitching  5°  and  frequency 
parameter  of  0.1,  it  was  something  which  could  be  visualised  from  a  steady  flow  pattern 
at  a  mean  angle  of  attack  0°.  The  quasi-steady  pressure  distribution  (Fig  5-41  and 
unsteady  pressure  distributions  (Figs  5-5  and  5-6)  were  lower  in  absolute  values  than  in 
any  of  the  above  cases,  and  they  formed  flat  curves  and  the  phase  lag  was  small. 
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5 . J  Comparisons  of  the  results  of  the  present  method  with  experiments  and  other  methods 

of  calculation 

(.Case  6)  A  flow  fit-id  having  an  extensive  supersoniv  lield,  terminated  by  a 
relatively  strong  shock  wave 

This  flow  field,  which  was  ol  the  same  type  as  that  used  in  case  2,  was  used  in  the 
calculations  for  the  NLR  7101  airfoil*  in  order  to  facilitate  comparison  with  the  experi¬ 
ments  (.M  =  0.7)  of  Tijdeman*.  Tijdeman  carried  out  experiments  witli  steady  flows  in 

relation  to  amplitudes  >  =  3.3°,  1.0°  and  1.5°  of  pitching,  and  obtained  the  pressure 

distributions  shown  in  Figs  6-4a  and  6-5a.  He  then  carried  out  experiments  with  unsteady- 
flow  about  an  airfoil  subjected  to  oscillation  in  pitch,  i(t)  =  3°  +  0.3°  sin  0.1841  , 
around  a  central  axis  IX  =  0.4,  Y  =  -1/60).  The  method  of  calculation  we  used  was  for 

C  0 

inviscid  flow,  but  the  experimental  results  will  have  the  influence  of  viscosity  and  the 
effect  of  the  wind  tunnel.  Accordingly,  steady  flow  calculations  were  made  for  M  =  0.7 
and  varying  angles  ol  attack,  the  angle  of  attack  a  =  l.5c  (Figs  b-4b,  6-5b  and  6-J)  , 
was  chosen  as  having  pressure  distribution  closest  to  the  results  of  experiments  carried 
out  at  a  =  3.0  .  The  value  1.3°  was  used  as  a  mean  angle  of  attack  in  the  calculations 
of  the  unsteady  flow  (Fig  6-3).  The  amplitude  of  displacement  of  C  during  one  period 
of  this  unsteady  flow  was  smaller  than  that  in  the  quasi-steady  flow,  as  will  be  seen  from 
Fig  6-4.  Quasi-steadv  pressure  distribution  and  the  unstable  pressure  distribution,  have 
a  high  pressure  peak  (Figs  6-5,  6-6,  6-7  and  6-8)  caused  by  shock  wave,  as  in  case  2. 

This  is  common  to  the  calculated  results  and  the  experiments.  However  this  peak  cannot  be 
predicted  from  the  solution  based  on  thin  airfoil  theory**,  also  shown  in  these  figures. 
The  results  of  calculations  made  by  the  present  method  are  in  good  qualitative  agreement 
with  experimental  results,  except  for  a  slight  deviation  in  the  position  of  the  shock 
wave  and  differences  in  absolute  values  of  steady  and  unsteady  pressure  coefficients. 

(Case  7)  Unsteady  phenomena  of  an  airfoil  under  design  conditions  of  being 
shock  wave-free 

The  supercritical  shock-tree  airfoil  NLR  7301  for  =  0.721  and  =  0,59^ 

has  been  designed  using  hodograph  theory.  In  the  experiments  made  by  Tijdeman,  in  which 

the  effects  of  viscosity  and  tunnel  walls  were  included,  values  M  =  0.744  and 

a  =  0.85  were  found  to  give  a  C  curve  approximating  the  C  curve  under  design  con- 

ditions.  In  the  experiment,  an  increase  in  angle  of  attack  by  0.5  produced  a  shock  wave 

on  the  upper  surface,  and  a  decrease  by  0.5°  produced  no  shock  wave  but  two  peaks  in  the 

C  curve.  As  will  be  seen  from  the  above,  a  slight  change  in  angle  of  attack  produced  a 
P 

great  change  in  steady  pressure  distribution  on  the  upper  surface  but  only  a  small  change 


*  This  airfoil  had  been  designed  by  the  hodograph  method'  of  Boerstoel  for  obtaining  a 
flow  free  of  shock  waves.  Being  relatively  thick  (maximum  thickness  16.5%)  and  blunt- 
nosed,  this  is  considered  to  be  representative  of  modern  supercritical  airfoil  designs 
Because  of  its  shape,  the  hypothesis  of  small  disturbance  is  liable  to  be  violated  and 
a  small  disturbance  equation  cannot  be  useu.  The  full  potential  equation  or  the  Euler 
equations  should  be  used.  The  Cartesian  grid  system  proposed  by  Carlson  and  used  by 
Isogai7  is  unsuitable  for  this  type  of  airfoil  with  a  blunt  nose. 

**  It  has  been  shown  by  experiments  that  an  equation  approximating  a  thin  airfoil  can  be 
used  satisfactorily  for  making  predictions  relating  to  subsonic  and  supersonic  flows, 
but  it  is  known^  that  the  equation  does  not  closely  agree  with  experimental  results  in 
the  case  of  supercritical  flow. 
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on  the  lower  surface  (.Fig  7-4a) .  In  these  calculations,  values  of  M _  =  0.7JI  and 

<  =  -0.1°°  were  found  to  give  an  approximation  to  the  C  curve  under  the  design  condi- 

o  ^ 

tuns.  This  resultant  angle  of  attack  a  =  -0.19  was  used  as  the  standard  from  winch 
steady  pressure  distributions  were  obtained  when  the  angle  of  attack  was  increased  or 
decreased  by  0.5°  (Fig  7-Ab)  .  As  is  seen  from  the  isobaric  contours  shown  in  Fig  7-J, 
a  small  change  in  the  angle  of  attack  results  in  a  large  change  in  the  flow  field  from 
tlte  design  conditions.  The  quasi-steadv  pressure  distribution  (Fig  7-5)  ,  has  a  broad 
bulge  on  the  upper  surface  in  both  experimental  results  and  the  results  of  the  present 
calculation.  The  curve  from  thin  airtoil  theory  differs  qualitatively  from  these, 
indicating  that  thin  airfoil  theory  is  unsuitable  for  prediction. 

Unsteady  phenomena  resulting  from  the  choice  of  i  =  0.5°  and  k  =  0.181  and 
the  design  conditions  are  as  follows.  The  unsteady  pressure  distribution  (Fig  7-3)  on 
the  lower  surface  became  a  relatively  smooth  curve  and  there  was  good  agreement  between 
the  experimental  results  and  the  results  of  the  present  calculation.  On  the  upper  surface, 
the  shape  was  similar  but  there  was  a  quantitative  difference.  Referring  to  Fig  7-7, 
curves  for  absolute  values  of  the  unsteady  pressure  distribution  coefficient  obtained  in 
the  experiment  and  the  present  calculation  have  a  peak  resulting  from  periodic  formation 
(Fig  7-3)  of  a  weak  shock  wave  and  the  phase  curves  are  smooth  for  about  40^50%  of  the 
chord  length  and  then  they  turn  upwards  sharply  owing  to  the  presence  of  the  shock  wave. 

As  in  the  case  of  quasi-steady  flow,  the  results  of  thin  airfoil  theory  differed  con¬ 
siderably  from  those  of  the  experiment  and  the  present  calculation.  As  is  seen  from  the 
above,  the  results  of  using  our  method  based  on  the  inviscid  theory  differ  quantitatively 
from  those  of  experiment,  but  qualitatively  they  capture  principal  characteristics  of  an 
unsteady  flow  about  an  oscillation  about  the  shock-free  design  condition. 

Ssruvk:  The  wind  tunnel  used  by  Tijderaan  in  his  experiment  was  capable  of  producing 

a  relatively  low  Reynolds  number  of  only  3.J  '  10fe.  This  signifies  that  the  transition 
from  laminar  flow  to  turbulent  flow  in  that  wind  Lunnel  takes  place  downstream  oi  that  in 
free  flight  (Re  =  30  '  50  v  10*).  Tijderaan  also  carried  out  an  experiment  with  a  transi¬ 
tion  zone  fixed  at  30%  chord,  which  is  similar  to  that  in  free  flight.  Under  the  condi¬ 
tion  of  a  strong  shock  wave  in  case  6,  there  was  no  effect  of  a  transition  zone  because 
spontaneous  transition  occurred  at  about  the  same  position,  but  it  will  be  seen  from 
Figs  7-8  and  7-9  that  a  flow  was  very  sensitive  to  the  transition  zone  under  the  design 
condition  in  that  case.  It  will  also  be  seen  from  the  figures  that  the  results  of  the 
present  calculations,  which  did  not  include  the  effects  of  viscosity  or  turbulent  flow, 
were  closer  to  the  results  of  experiments  with  a  transition  zone  than  to  those  of  experi¬ 
ments  in  spontaneous  transition  without  a  transition  zone.  It  is  difficult  to  say  that 
the  flow  field  in  this  experiment,  undesirable  phenomena  of  the  formation  of  an  oblique 
compression  wave  from  the  transition  zone,  and  its  transformation  into  a  compression  wave 
on  being  reflected  by  a  sonic  line  were  observed.  All  the  same,  a  flow  under  the  design 
conditions  reacts  very  sensitively  to  a  slight  change  in  condition. 


(Case  8)  Unsteady  transonic  wave  accompanied  by  a  weak  shock  wave 

Davis,  carried  out  experiments  will)  the  symmetrical  airfoil  NACA  h4A0lU,  wi:l. 

a  fine  and  rounded  leading  edge,  with  free  stream  Mach  number  0.8.  Numerical  solutions 
for  it  were  sought  by  Magnus-’’"’,  using  the  Euler  equations,  on  assuming  inviscid  flow, 
and  by  Isogai"1  who  used  a  full  potential  equation  and  a  transonic  small  disturbance  <TSDi 
equation.  Similarly,  we  sought  three  steady  (lows  under  the  conditions  of  the  same  free 
stream  Mach  number  and  angle  ol  attack  ot  0°  •  1°.  Also  we  obtained  unsteady  flows  at  the 
mean  angle  of  attack  0°  and  in  amplitude  of  pitching  I ^ ,  for  comparison.  Isobaric  contours 
and  C  curves  tor  steady  and  unsteady  flows  obtained  by  the  present  method  of  calculation 
are  shown  in  Figs  8-2,  8- i  and  8-4.  The  airfoil  being  symmetrical  and  having  a  mean  angle 
of  attack  ot  0°,  the  flow  field  at  a  certain  cycle  in  Fig  8-1  is  the  same  as  the  inversion 
of  a  flow  field  at  a  cyele  differing  by  0.5  in  absolute  value.  The  steady  flow  shown  in 
Fig  8-2  (Fig  3--<)  has  no  supersonic  region  (shook  wave  absent)  on  the  lower  surface  when 
the  angle  of  attack  is  l1,  but  the  unsteady  flow  shown  in  Fig  8-2  (Fig  8-4)  has  supersonic 
regions  (a  shock  wave  present)  on  both  upper  and  lower  surfaces  at  all  cycles.  t'ompari- 
sons  of  the  results  of  the  experiment,  the  present  calculation,  and  other  calculations 
are  shown  in  Figs  8-5  to  8-8.  In  the  steady  pressure  distribution  (Fig  8-5)  at  angle  of 
attack  0°,  the  shock  wave  according  to  the  present  calculation  is  located  upstream  of  the 
others,  but  its  width  is  the  smallest,  followed  by  a  fall.  The  present  method  of  calcula¬ 
tion  captured  the  shock  wave  sharply.  The  peak  resulting  from  the  periodic  displacement 
of  the  shock  wave  i n  the  unsteady  pressure  distribution  were  obtained  in  decreasing  order 
of  size  (Figs  8-7  and  8-81  from  Magnus'  method  using  the  Euler  equations,  the  present 
method  using  the  same  equations,  Isogai's  method  using  a  full  potential  equation,  and 
Isogai 's  method  using  a  small  disturbance  potential  equation.  Considering  the  characteris¬ 
tics  o‘‘  these  equations,  this  order  seems  appropriate.  Obtaining  different  results 
(Fig  -6)  from  using  the  same  KuK  r  equations  resulted  from  differences  in  grid  system  and 
the  finite  difference  method.  Among  the  methods  used,  those  which  were  particularly  close 
to  unsteady  experimental  values  were  the  present  method  and  Isogai's  method  using  the  lull 
potential  equation.  The  results  from  both  methods  approached  the  experimental  values, 
thus  they  both  seem  to  be  satisfactory  for  prediction  in  this  case.  This  case,  in  which 
the  airfoil  thickness  is  small  and  the  leading  edge  is  sharp,  is  also  particularly 
suitable  for  Isogai's  method  using  the  small  disturbance  potential  equation  in  which  the 
airfoil  slope  is  handled  on  a  flat  plate  or  the  method  of  Isogai  using  a  full  potential 
equation  employing  a  Cartesian  grid  system  which  is  considered  to  be  weak  when  used  for 
thick  and  rounded  leading  edges. 

'"Kre:  Because  of  the  use  of  a  maximum  number  of  grid  points  of  121  v  11  in  the 

program,  no  satisfactory  solution  could  be  expected  in  this  case.  Accordingly,  a 
A\/4-wide  grid  prepared  for  use  only  in  the  vicinity  of  vertically  usyrame t r iea 1  shock 
wave  was  used  vertically  symmetrically  in  the  vicinity  of  the  leading  edge. 


i  Required  computing  time 


The  computing  times  required  tor  the  calculation  ot  numerical  values  in  sections  b.l 
and  l.d,  using  a  FA  COM  J30-75  computer  |C1T  (ordinary  computer)  +  Al’l'  (.array  processor)!, 
in  accordance  with  the  routine  described  in  section  A  are  as  follows: 

(I)  the  formation  of  a  single  grid  and  the  about  db  seconds  with  C1T 

computation  of  associated  differential 

coef f icient s 

(J) ,l3)  computation  of  one  steady  flow  about  30 "  90  seconds  with  AIT 

(-»)  the  processing  of  results  pertaining  to  a  about  9  seconds  with  ('IT 

quasi-steadv  flow 

lb)  computation  of  an  unsteady  flow  cycle  about  d5  "  J  seconds  with  AIT 

lb)  processing  of  results  pertaining  to  an  about  dO  seconds  with  HIT. 

unsteady  flow 

Most  of  the  computing  times  were  used  in  performing  routines  Id)  and  ( 3) ,  using 
program  FLOW.  In  the  computation  of  steady  flows  under  (d)  and  (1),  the  time  requirement 
varied  greatly  according  to  the  division  of  time  chosen  for  the  selection  of  initial 
values  and  for  maintaining  the  complexity  and  stability  of  flow  patterns.  I'h  processing 
of  results  pertaining  to  an  unsteady  flow  in  (.3)  is  similar,  but  it  depends  on  :  requester 
parameter  and  is  independent  of  the  selection  of  initial  values.  The  average  total  time 
required  in  the  case  of  computing  for  quasi-steady  flow  and  four  cycles  ot  unsteady  slew 
was  about  5  hours. 

s . w  Deviation  from  Rankine-Hugoniot  relation 

The  following  Rankine-Hugoniot  relation  holds  theoretically  in  a  steady  normal 
shock  wave  in  an  inviscid  flow 


—  =  I  +  (>r  -  l)  . 

b  +  l  V  i  / 

In  the  above  equation,  the  subscripts  1  and  d  indicate  positions  immediately  before  and 
alter  the  shock  wave.  The  suitability  of  a  particular  numerical  solution  as  a  solution 
for  an  inviscid  flow  can  be  judged  by  the  degree  ol  its  closeness  to  the  Rankine-Hugoniot 
relation. 

With  a  view  to  finding  the  suitability  of  the  present  method  of  calculation,  the 
quasi-steady  and  steady  flows  (13  in  number)  which  could  be  considered  to  contain  normal 
shock  waves  were  selected  from  among  those  obtained  in  the  examples  of  numerical  calcula¬ 
tions  referred  to  earlier.  Pressure  rise  at  a  shock  wave  was  obtained  from  Mach  number 

and  pressure  obtained  at  the  actual  mesh  points.  These  values  are  plotted  together  with 

’3 

results  of  other  methods  of  calculation  in  a  diagram  (Fig  9)  devised  by  Lomax,  ,  ‘  . 

As  expected,  the  solutions  by  time-dependent  methods  based  on  the  Euler  equations,  includ¬ 
ing  the  solution  by  the  present  method  of  calculation  were  found  to  better  satisfy  the 
Rankine-Hugoniot  relation  than  solution  by  a  relaxation  method  based  on  a  potential 
equation. 


However,  bi’i  ausi'  ol  the  presence  ol  a  boundary  layer  in  the  experiment,  the  strong; 
>'l  the  shock  wave  measured  on  an  airfoil  is  generally  lower  than  what  can  be  predicted 
t rom  the  above  relation.  Because  of  this,  a  numerical  solution  using  a  potential  equal ie 
which  is  Lower  in  degree  of  approximat ion  than  Euler  equations  for  an  inviscid  flow  car. 
often  be  closer  to  an  experimental  result.  For  obtaining  a  solution  justifiably  closer 
to  an  experimental  value,  it  is  necessary  to  solve  Navier-Stokes 1  equation  tor  viscous 
f  low . 

h  CONCLTSION 

This  work  proposes  a  method  for  the  calculation  of  an  inviscid  compressible  flow 
about  a  :  wc-d  imens  i  o;ia  1  airfoil  oscillating  sinusoidally  in  pitch.  With  a  view  to  hata.li 
accurately  at:  airfoil  even  with  a  rounded  and  thick  leading  edge  (such  as  the  NLR  Viol 
airfoil),  which  appears  difficult  to  solve  by  previously-published  methods  for  :alcu 1  at  in 
unsteady  flows.  A  coordinate  transformation  was  devised  and  a  computational  plant  formed 
The  Euler  equations  which  have  no  approximation  were  used  to  model  an  inviscid  flow,  and 
the  iinice  difference  method  of  Lax-Wendroff  tcith  artificial  viscosity  was  used  for  its 
solution.  Then,  the  effectiveness  of  the  present  method  was  examined  by  carrying  out 
numerical  ca lcul at  ions ,  using  four  airfoils  (Joukowski,  NACA  0012,  MLR  730!  and 
NACA  n 4 AO  10  airfoils)  within  a  wide  range  (C.7  '  1.4)  of  Mach  number.  Consequently,  we 
succeeded  in  numerical lv  capturing  various  types  of  unsteady  shock  wave  phenomena 
(processes  of  generation,  growth  and  disappearance  of  a  weak  shock  wave,  processes  of 
propagation  of  a  strong  shock  wave,  a  'fish  tail'  shock  wave  system  and  its  process  of 
decaying,  the  process  of  propagation  of  a  shock  wave  separated  from  the  leading  edge,  and 
others).  During  this  time,  we  proved,  as  was  expected,  that  the  discontinuities  across 
captured  shock  waves  satisfy  better  the  Rankine-Hugoniot  relations  than  do  those 
calculated  using  the  fuel  potential  equation.  The  present  method  of  calculation  seems 
fully  workable  for  predicting  inviscid  and  unsteady  aerodynamic  phenomena  but  is 
unsuitable  for  unsteady  flow  with  a  strong  viscous  interference.  If  a  more  powerful 
computer  than  is  now  available  became  available  in  the  near  future  then  it  should  be 
possible  to  get  numerical  solutions  for  unsteady  flows  with  strong  viscous  interference 
by  modelling  Che  viscous  flow. 

The  author  is  grateful  to  Technical  Officer  Isogai  for  participating  in  the 
discussion  and  for  supplying  the  data. 
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Figs  1-3  -  1-6 


Fig  1-3  Steady  and  unsteady  pressure  coefficients 
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Fig  1-4  Quasi-steady  pressure 
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Fig  1-5  Unsteady  pressure 
distributions 
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Fig  1-6  Unsteady  pressure  distributions 


Figs  2-4 


Fig  2-4  Unsteady  isobaric  contours  k  =  0.25 


Fig  2-5  Steady  and  unsteady  pressure  coefficients 
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Fig  3-3  Unsteady  isobaric  contours 
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Fig  3-7  Unsteady  pressure 
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Fig  6-4  Steady  and  unsteady  pressure  coefficients 


QUASI  -STEOTY  FLOW 


Fig  6-5  Comparison  of  quasi-steady  pressure 
di stributions 


PRESENT  CALCULATION 


NLR  7301  UPPER  S URfftCE 
*m- 0.70  4«=0 .S' 

■  o  =  3.0*  k=0 . 192  flCp. 


—  EXPER  I n  E  N  T 1  T I  JOE  HAN  ) 
thin-airfoil  theory 


Fig  6-6  Experimental  distributions  of 
unsteady  pressure 


PRESENT 
•■calculation 
n„=o  .  7  0 
«o  =  l  -S’ 

*  A  «  =  0 . 5  ’ 

'c  K=0.i92 

-V 

•0 

■  m 
x 


I-OWe’r  SURFfiCE* 


I  .0 


NLR  7301 


Fig  6-7  Unsteady  pressure  distribution 
by  the  present  method  of 
calculation 


Fig  6-8  Comparison  of  unsteady  pressure 
distributions 


Fig  7-4  Steady  and  unsteady  pressure  coefficients 
ouasi-STEno'f  f'jJw  NCR  7301 


°  €xP£Rl*E*TlT!jOE"«Ni  0  PRESEN'  C  *  L  C  •.  -  ®  M  3  S 


Fig  7-5  Comparison  of  quasi-steady  pressure  distributions 


Fig  7-6  Comparison  of  unsteady  pressure  distributions 


(  q  )  l »p£« I ntHJ 
(  1 J  J0f  ) 
r  m- 0  •  7 46 
« o  :  0  ■  85  ‘ 
*■*  *  4«-0  5  ' 
a  h  ;  0  •  1 8 1 


(D  »  P«t  51 W  T 

■  CHI C  Ut  «  T  I  ON 

nmz o .72 i 

•  O  '-  -  0  .  i  9  ‘ 

\  i§zO-S‘ 

*  =  0-181 


T  H  |  N-  H  I  *f  o  f  L 
T  M  F  Q  R  » 


*  ♦  RESENT 

caxuiATiw 
n.=0-721 
„  i0=-<M9- 

*«=o. s* 

k=0.;8j 


NLR  7301  oo  oo  it 

Fig  7-7  Comparison  of  unsteady  pressure 
di stributicns 


Co*** 


exrm’OT  „  _0.74S 
(TlJOEf***)  .‘=0.8S- 
onwtsincw  strip  4»=o.s' 
awinjw*.  rrawsinoH  k=c-i8i 


* . 


C.  »D 
-o  *’  c 

<"  ’»  / 

*V 


txrwncm  «0=o  es‘  k=o.iei 

( rixowi 

mnm.  nwKstnoN 

OH.=fl-744  TRfWSITIW  STRIP 
4 4=0 -56'  «  n.=fl-742 

S,  •  N„=0.743  4«=0-48' 

44=0-76* 


Fig  7-8  Comparison  of  unsteady  pressure 
distributions 


/  j 

<«/  / 

f4i %y,.X 


a  t  ,/ 

- 

s-  v 

li  >* 

f  I  P«swr 

'  1  CfLOJLflTIOH 

.  Ofi-  « 
o  i  •,=-0.19' 

}  4«=0  -5* 

k  =0.161 


.  ♦ 
\.®*V 


Fig  7-9  Comparison  of  unsteady  pressure 
distributions 


[tillin' 


iiiilTr 


Figs  8-4  -  8-6 


Fig  8-4  Steady  and  unsteady  pressure  coefficients  obtained  by 
the  present  method  of  calculation 
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Fig  8-5  Comparison  of  steady  pressure 
coefficients 
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Fig  8-6  Comparison  of  quasi-steady 
pressure  distributions 
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Fig  8-7  Comparison  of  unsteady 
pressure  distributions 


Fig  8-8  Comparison  of  unsteady 
pressure  distributions 
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Fig  9  Deviation  from  the  Rankine-Hugomot  relation 


